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Abstract. We show carriers dynamics on bulk Silicon in the sub pico-second time scale. The 
results are obtained from a first-principles implementation of the the Kadanoff-Baym equations 
within the generalized Baym-Kadanoff ansatz and the complete collision approximation. 
The resulting scattering term is similar to the scattering described within the semi-classical 
Boltzmann equation pQ. 


1. Introduction 

The description of carriers dynamics in semi-conductors is crucial for a correct modelization of 
electronic devices. Silicon (Si) in particular is the fundamental building block of semi-conductors 
physics and microelectronics industry [2J. 

The reference equation used in the literature is the semi-classical Boltzmann equation (sc-BE) 
which describes the scattering of electrons and holes with lattice vibrations, electron-phonons 
(e-p) scattering, with other carriers, electron-electron (e-e) scattering and possibly with defects. 
Its implementation is usually based on models, where the characteristic of the system are taken 
into account with ad-hoc parameters. In the present work we want instead to describe such 
scattering events within a first-principles approach. 

Also the sc-BE assumes a semi-classical regime, as opposite to the quantum-kinetic regime 
which should, in principle, describe the electrons propagation. However the miniaturization of 
Si-based devices to the nano-scale regime and the never ending search for faster devices call for 
a deep understanding of the fundamental quantum-mechanical process that governs the carriers 
dynamics 13 i- Moreover the recent development of ultra-short laser pulses has triggered the 
necessity to understand the quantum kinetics of electrons and holes on the sub pico-second ( ps ) 
time-scale. 

The non-equilibrium Green’s functions (NEGF) approach offers a quantum-kinetic equation, 
the Kadanoff-Baym equation (KBE), which describes the quantum dynamics through the many- 
body Green’s function G < (t, t'). The solution of the full equation however is not feasible within a 


first-principles approach, and also for simple models, an approximate version is usually adopted. 
The most common approximation is the generalized Baym-Kadanoff ansatz (GBKA), which 
reduces the KBE to a closed equation for the time-diagonal Green’s function, i.e. the time- 
dependent density matrix of the system. If also the complete-collision approximation (CCA) 
is employed the scattering processes are described at the same level of the sc-BE jT], however 
with a smeared energy conservation due to possible violations on the ultra-short time regime. 
We will refer to this equation as the CCA-KBE. 

The CCA-KBE has been recently implemented in the yambo [5J code. Projected onto 
the Kohn-Sham (KS) wave-functions, it describes the evolution of a parameter free G-(£,£), 
which can be used to reconstruct the electronic occupations: fi = Here we used 

the generalized index i — {nkcr}. Although the scattering terms still have a semi-classical 
form, the CCA-KBE offers a suitable starting point for future development towards a real 
quantum kinetics equation. It also has other advantages. For example it can correctly model 
the interaction with ultra-short laser pulses. Indeed the first-principles CCA-KBE has been 
shown to reproduce the propagation of the polarization [ 6 ], P — J2i^j Gfj(t, t ) induced by a laser 
pulse. 

In the present work we explore the performance of its implementation on bulk Silicon and 
in particular we will focus our attention of the scattering term Si(t). The details and the 
approximations needed for a first-principles implementation of this term have been discussed in 
Ref. [7]. 


2. The theoretical method 

The theoretical approach used is described in Refs. urn- Here we report the main equations. 
The CCA-KBE equation reads 

- [Hmb + H RT (t), G<(t)]. j = Sij(t) ( 1 ) 

with 


Hdft [p eg ] = H 0 + V H [p eq ] + V xc [p eq ] ( 2 a) 

H M b[G = H DFT [p eq ] + (£[G<’ e<? ] - V xc [p eq \) (2b) 

H RT [G<(t)](t) = V H [Sp(t)\ + (E[<$G<(t)], G<’ eq }) + U ex \t) ( 2 c) 

5(t) = {E>,G<} + {E<,G'>} (2d) 


Hdft is the standard KS Hamiltonian. Hmb describes the quasi-particle (QP) corrections. 
Finally the Hfct Hamiltonian instead includes the external field and the variation in the Hartree 
potential and in the self-energy [ 6 J. In the present work we replace its role with an ad-hoc 
initial guess for the electrons and holes populations. 

The term S(t) describes the carriers scattering [7]. Since the electronic occupations are the 
diagonal terms of the G < ij 1 we focus our attention on the diagonal of eq. [l] The latter has 
only contributions from the S(t) term, since we replace the role of H^t with an initial guess 
and Hmb gives no contribution on the diagonal. Indeed the CCA-KBE, together with a QP 
approximation for the retarded and advanced Green’s function [7], has the form 


a./. (*)U. = —7, <e, («)/f(*) + 7’(*), 


(3) 


with fj e ^ = fi and = 1 — fi. The electron and hole lifetimes, and 7 ^, include both 

(e-p) and (e-e) contribution: ^ e / h ) = 7 ^ 6// ^| e _ e + 7^ h ^\ e _ p - 7 ^ 6// ^ = 7 ^ e ^| e _ e + 7^ h ^\ e _ p - 
In the e-e case we have that: 


J e / h ) 

H 


oc 5] \Wi^ kl \ 2 P (Ay - A kl ) if /e) (t) fi h/e) C t ) ft e/h) C t ) , 


(4) 


e—e 



1=0.0000000 fs 


1=10.500000 fs 


1=31.500000 fs 


1=63.000000 fs 










W 2961 1 K ^ 


5 -4 -3 -2 -1 0 1 2 3 4 5 

[eV] 

1=252.00000 fs 

V 

rW 4176 - 1 K 

l^-8137.7 K 

V 

5 -4 -3 -2 -1 0 1 2 3 4 5 

[eV] 

1=829.50006 fs 

n 

V T elecs= 31263 K 

\ 


-5 -4 -3 -2 -10 1 2 3 4 5 

[eV] 


Figure 1 . Carriers occupations distribution as a function of time 


while in the e-p case 


3 A 


(5) 


In Eqs. [4][5] Aij — — ej and A represents a generic phonon mode with momentum q and branch 

g. N+ = (1 + N\ (T) /2) and N* — N\ (T) /2, with N\ ( T ) the Bose distribution function at 
energy uo\ and temperature T. In the present work we assume that the phonons distribution 
does not change during the dynamics. W is the statically screened component of the Couloumb 
interaction whose frequency dependence is embodied in the P functions. These represent a 
smeared energy conservation condition. Finally, the g are the matrix elements of the screened 
ionic potential derivative, calculated within DFPT [8, 9J. 


3. Results on bulk Silicon 

In the present work we used the local-density approximation (LDA) to describe bulk Silicon 
with norm-conserving pseudo-potentials with an energy cet-off of 20 Ha. For the Brilloine 
zone we use a uniform 4x4x4 sampling and we consider the first 7 bands, i.e. 4 occupied bands 
and 3 empty bands. The QP corrections are then described with a scissor operator of 0.8 eF, 
which is known to be a very good approximation for buk Si. Then we compute the phonons and 
the e-p matrix elements with a DFPT run using the same parameters. The CCA-KBE is then 
propagated in the KS wave-functions space, but on an auxiliary 15x15x15 k-grid. This bigger 
grid is required to ensure a proper description of the e-p scattering. The KS eigen-values are 
computed on this finer grid, while all other quantities are interpolated from the 4x4x4 grid with 
a constant nearest-neighbor interpolation. The phonons are assumed to remain in equilibrium 
according to a Bose distribution of 300 K. The time step for the propagation is dt = 10 as, 
while we use a second-order Runge-Kutta integration scheme. The conservation of the number 
of carriers during the dynamics is guaranteed by the detailed balance, as discussed in Ref. [7]. 






































We consider a starting point which mimics a high number of electrons excited from a laser 
pulse centered at the optical band-gap of Si. Thus we create near the V point, a packet of holes 
in valence and a packet of electrons in conduction. We use here a very high density of carriers, 
which would physically lead to a drastic change to the characteristics of the material. The choice 
is done to test the implementation in a high-carriers regime. Indeed here we are interested in 
exploring two selected physical effects, i.e. the e-e scattering and the e-p scattering, in this 
regime. In a lower density regime the effect of the e-e scattering would become smaller and it 
would be very difficult to test the reliability of the implementation and eventual errors. 

In doing this we will neglect some physical effects, such as the softening of the phonon 
frequencies, the change in the band structure due to the extra carriers, the deviation of the 
phonons distribution from a Bose function and the sub-sequent heating of the lattice. All these 
effects becomes negligible in a lower density regime whereas instead the weaker e-e scattering 
would still have a role. 

In Fig. [l] we see different time-snapshots of the carriers distribution as a function of the 
energy. We see that the electrons and the holes, which at the initial time are in a highly non¬ 
equilibrium situation, relax towards two Fermi distributions with changing temperatures. In 
the holes channel, carriers are closer to their minimum (i.e. the valence band maximum, VBM) 
and they reach a lower temperature Fermi distribution. In the electrons channel, instead, we 
can appreciate the different role played by the e-e scattering and the e-p scattering. Phonons 
have, on average, an energy of few meV and can scatter electrons only to levels energetically 
nearby, the result is a packet of electrons which keeps approximately the initial spread in energy 
and slowly moves towards the conduction band minimun (CBM). This is poorly visible because 
meanwhile the e-e channel is able to directly scatter electrons to the CBM. 

The resulting temperatures, obtained as 
the fit of the two Fermi distributions, 
are reported in Fig. [2| Here, together 
with the electron and hole temperature 
for the simulation with both e-e and e-p 
scattering on we report also the fit in case 
only one of the two relaxation channels is 
open. The e-e scattering (red lines) makes 
the electrons and holes interact and thus 
increase the holes temperature lowering the 
electrons temperature. However once the 
two sub-systems reach a thermal equilibrium 
nothing more can happen. Instead the 
e-p channel (blue lines) makes the carriers 
exchange energy with the lattice. Thus 
both electrons and holes decrease their 
temperature. However since phonons do not 
have enough energy to cross the electronic gap, the two processes are independent. Finally, 
when both scattering channels are open the resulting effect is a combination of the two. In this 
high carriers density regime, the e-e scattering dominates the dynamics in the initial part of the 
dynamics while, once electrons and holes reach a similar temperature, the e-p channel becomes 
the dominant channel. 

The temperatures reached by electrons and holes depends mainly on the initial guess (or in 
case a laser pulse were used on its characteristics) and on the Silicon band structure. They are 
instead almost independent on the total number of carriers created. Indeed lowering the carriers 
density the main change is that the e-e scattering become less effective while all occupations 
numbers are rescaled. In the limit of very low carriers concentration the evolution of the carriers 



Figure 2. Overview of the electrons and 
holes temperatures for different scattering 
mechanisms 
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Figure 3. Carriers occupations on the band structure as a function of time 


temperature becomes very similar to the case with only the e-p channel active (data not shown). 

Finally in Fig. [3] the carriers dynamics is represented on the whole band structure. Here we 
see from another perspective the difference between the two channels. Holes remain basically 
close the VBM, just cooling down, while electrons move along the band structure going from T 
to the minimum, nearby X. 

In conclusion we have presented a first-principles implementation of carriers dynamics in 
bulk Silicon. The most critical aspect of the present simulations is the use of a double grid 
technique for the sampling of the Brilloine zone to guarantee a correct description of the energy 
transfered from electrons to phonons. The global picture reveal a three systems dynamics: the 
electrons in the conduction band, the holes in the valence band and the lattice. The three 
systems exchange energy, electrons an holes vie the e-e scattering, while carriers and phonons 
via the e-p scattering. Both electrons and holes reach an internal meta-stable state, i.e. they 
are distributed according to a Fermi distribution, very quickly (approx 0.1 ps). However the 
complete equilibration, either in between electrons and holes (less then 1 ps) or with the lattice 
(not yet reached after 1.5 ps), takes longer time. This is, to the best of our knowledge, the first 
fully ab-initio and parameters free description of carriers dynamics in Silicon. 
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